Ion-assisted ground-state cooling of a trapped polar molecule 
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We propose and analyze a scheme for sympathetic cooling of the translational motion of polar 
molecules in an optical lattice, interacting one by one with laser-cooled ions in a radio-frequency 
trap. The energy gap between the excitation spectra of the particles in their respective trapping 
potentials is bridged by means of a parametric resonance, provided by the additional modulation of 
the RF field. We analyze two scenarios: simultaneous laser cooling and energy exchange between the 
ion and the molecule, and a scheme when these two processes take place separately. We calculate 
the lowest final energy of the molecule and the cooling rate depending on the amplitude of the 
parametric modulation. For small parametric modulation, the dynamics can be solved analytically 
within the rotating wave approximation. 
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I. INTRODUCTION 

Ultracold molecules are the focus of increasing atten- 
tion both in physics and in chemistry [IHH- On one 
hand, recent experiments have succeeded in creating cold 
molecules via techniques like photoassociation of ultra- 
cold atoms [i| or magnetically tunable Feshbach reso- 
nances and in cooling them to the ground vibra- 
tional state, e.g. via optical pumping [^-Q. On the other 
hand, theoretical proposals have shown the relevance of 
such systems, especially heteronuclear molecules bearing 
nonzero electric dipole moment, for a range of applica- 
tions from metrology [ToL [Tl| to quantum information 
processing [HI, 03 and the simulation of quantum many- 
body systems |14| - {l61 |. 

The preparation of trapped polar molecules in the 
ground state of both their internal rovibrational and 
center-of-mass motional states is therefore of crucial im- 
portance for further scientific and technological develop- 
ments. Different approaches have been put forward, from 
the coupling to cavities, in the optical [13, 03 as well 
as in the microwave regime [l^] to laser cooling |20l . l2l| 
and sympathetic cooling via collisions with neutral atoms 
[iH [23| • Cavity- based approaches, despite significant ex- 
perimental requirements, hold the promise of reaching 
the ground state also for the molecular center-of-mass 
motion. On the other hand, buffer-gas cooling schemes 
based on neutral "coolant" atoms are so far limited in 
the reachable temperatures - in the case of He [23| , sim- 
ply by the buffer gas' own temperature, and in the case 
of Rb 22] by losses through inelastic collisions. Typi- 
cally, different cooling schemes are applied to different 
temperature regimes of relevance for different types of 
experiments. The method we propose here is suitable 
for application in the ultracold regime in order to reach 



and maintain the molecular vibrational ground state via 
a sympathetic cooling procedure. 

Our approach builds on recent advances in the com- 
bined manipulation of charged and neutral particles in 
AMO physics, where e.g. cold collisions between ions and 
neutral atoms have been observed [24j and trapping of a 
single-ion within a Bose-Einstein condensate has been re- 
ported, including sympathetic cooling of the ion motion 
[25{. The scenario we envision takes one step forward 
from these pioneering experiments and involves a one- 
dimensional array of molecules trapped in an optical lat- 
tice and interacting with a chain of ions in a nearby ion 
trap. We assume the molecules to be in the deep Mott- 
insulator phase, so that tunneling between neighboring 
lattice sites is strongly suppressed. In this regime, the 
problem can be reduced to considering pairwise interac- 
tions between one molecule and one ion at a time along 
the direction transverse to the trap axis. 

In the following we will therefore focus on a single 
ion-molecule pair, and study in detail their interaction, 
with particular attention to their motional dynamics. In- 
deed, our idea for cooling the center-of-mass motion of a 
trapped molecule is based on a novel concept: the trans- 
fer of vibrational excitations to a neighboring, indepen- 
dently trapped ion that is subject to laser cooling. Re- 
peated application of this procedure allows the molecule 
to reach the trap ground state over a time scale compara- 
ble to typical ion cooling rates. Moreover, the procedure 
is entirely insensitive to the molecule's internal state, and 
therefore is suitable for applications where the latter en- 
codes information whose coherence needs to be preserved 
throughout motional cooling. The same applies, for in- 
stance, to atoms in an optical lattice used for quantum 
information processing, where quantum operations leave 
behind some residual heating of their external degrees 
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of freedom. The ion-atom interaction shows the same 
distance dependence as the ion-molecule potential dealt 
with in this paper, and therefore our proposal is appli- 
cable also to the latter type of scenario, allowing to use 
one or more trapped ions to sympathetically cool the mo- 
tion of qubit-carrying neutral atoms. This is relevant to 
ongoing experiments on ion-atom interactions that are 
being setup and carried out at several groups worldwide 

Mm- 

The paper is structured as follows. In Sect. [TT] we 
introduce the setup we have in mind and its descrip- 
tion in terms of two harmonic oscillators coupled by the 
long-range ion-atom interaction, modeled according to a 
quantum-defect approach. In Sect. Mil we present the 
simplest scheme of alternating cooling and energy trans- 
fer stages by exploiting a parametric resonance obtained 
with trap frequency modulation. In Sect. |lV]we discuss 
a different scheme in which the transfer and cooling pro- 
cesses take place simultaneously, so that a global master 
equation needs to be derived, which we do both within 
and beyond the rotating wave approximation in order 
to deal with arbitrary amplitudes of the parametric trap 
modulation. Finally, Sect. fVl contains a summary of our 
results and of future outlooks. 
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FIG. 1. Specific considered setup consisting of a trapped po- 
lar molecule and an ion interacting via long-range polarization 
forces. 

where the label i (to) refers to the ion (molecule) , p and r 
are the momentum and position operators, d denotes the 
distance between the traps, and to to , w m (to^, Ui) are the 
molecule (ion) mass and trapping frequency, respectively. 
Without loosing generality we can assume that d = de z , 
and for simplicity we have assumed spherically symmetric 
trapping potentials. 

B. Modeling as two coupled oscillators 



II. THE MODEL 



A. Setup 



As discussed above, we reduce the problem to con- 
sidering a single spinless polar molecule in its electronic 
and rovibrational ground state, with zero projection of 
the total angular momentum on the internuclear axis 
(A 1 E,^ = 0, J = 0). The molecule is confined in a 
far-off-resonance optical trap, created in a laser standing 
wave (an optical lattice), and it interacts with a single 
ion confined in a radio-frequency (RF) trap (see Fig.Q]). 
Because of the different nature of the trapping potentials 
for the two particles, we can assume that the distance 
d can be arbitrarily controlled by adjusting the external 
fields. A microscopic analysis of the ion-molecule sys- 
tem, presented in Appendix [5] confirms that for typical 
parameters of the radio-frequency and far-off-resonance 
optical traps one can neglect the cross effects of the 
RF and optical fields on the molecule and ion, respec- 
tively. Moreover, the derivation shows that at large dis- 
tances the ion attracts the polar molecule with a poten- 
tial V(r) -C 4 /r 4 , (r ->■ oo). Here, C 4 = d 2 m e 2 /(<oB), 
where d m is the molecule permanent dipole moment, e 
is the ion charge, and B denotes the molecule rotational 
constant. An effective Hamiltonian for the system can 
be written as 



2 m, . 



1 2 2 

-m v u) v T v 



V(\r l -r m +d\), 



We assume that the distance d between the two traps is 
sufficiently large, such that the effects of the ion-molecule 
interaction can be described perturbatively as distortions 
of the trapping potentials. In this regime the molecule 
and the ion remain well separated in the course of dynam- 
ics, and the ion-molecule interaction can be described by 
the asymptotic part of the potential: V(r) — —C^/r 4 . 
This allows also to neglect the effects of the tunneling 
through the potential barrier separating both particles, 
thus avoiding inelastic collisions. An expansion of V(r) 
up to second order terms in the distance r leads to 
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+ 20-jg(zi - Zi)(z m - z m ) - 4—gp m pi, (2) 

where u>l = uj 2 . — 20C 4: /(m l ,d e ) and ib 2 — uj 2 + 
4C4/ '(m,,d 6 ) are the trapping frequencies modified by the 
ion-molecule interaction, z v denotes the equilibrium po- 
sition of the particles and p = (x, y, 0). The Hamiltonian 
([2) describes a system of two coupled harmonic oscilla- 
tors. Since the motion in the spatial directions x, y and 
z separates, in the following we focus on the dynamics 
along z, noting that a similar analysis can be carried out 
for the x and y directions. Introducing the usual annihi- 
lation and creation operators: 
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with Qi = Zi - Zi, q m = z m - z m , pi = (pi) z and 
Vm = (Pm)z> we rewrite the part of the Hamiltonian @ 
describing the motion in the z direction in the following 
way: 

H z = hu> m a'a + huiib^b + huj c (a + ar)(b + b'), (5) 

where we omit the ground-state energy term. Here, u> c 
denotes the coupling frequency: ftiu c — lOCilmk/d 6 , and 
l v = -iJh/Jm^Qv), v — i.m are the harmonic oscilla- 
tor length for the ion and the molecule respectively. The 
value of oj c is bounded by the condition that the quadratic 
terms in the expansion of the potential are positive, 
which yields the stable equilibrium position for both par- 
ticles. This is fulfilled for C \d~^ '(1 / 'mauf + l/m m uj^ n ) < 
(|) 6 /20 ~ 0.017, which brings some maximum value of 
the coupling frequency 
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For typical parameters for polar molecules in optical po- 
tentials and ions in RF traps, w™ ax is of the order of 
10kHz. Considering an ultracold KRb molecule in the 
rovibrational sing let ground-state X 1 ^ (d m = 0.566D 
and B = 1. 1139GHz Q) confined in an optical trap of 
frequency tu m = 2ir x 100kHz, and a 40 Ca + ion confined 
in an RF trap of frequency = 2tt x 1MHz, we ob- 
tain cj" lax = 27.3kHz corresponding to a trap separation 
d = 327nm. 

The Hamiltonian ([5]) can be diagonalized by apply- 
ing the transformation described in Appendix [Bj which 
yields 



H z = hajia'a + hu}20 /3, 



(7) 



where u)% and u>2 are the frequencies of the two eigen- 
modes 
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with the upper (lower) sign referring to the frequency 
u)\ (W2). For a weak coupling between the oscillators 
(to c -C \u!i — uj m \), the mode cva (/3'/3) is localized mainly 
on the molecule (ion). At distances when the coupled 
oscillator model remains valid, u), ~ LOi and <D m ~ uj m , 
and throughout the rest of the paper, we will neglect 
small corrections to the trapping frequencies due to the 
interaction potential, using the bare trapping frequencies 
uii and u)m . 



C. Quantum-defect treatment 

The validity of the coupled-oscillators model ([5]) can 
be verified by comparing it with a calculation based 
on a quantum-defect modeling of the short-range part 
of the interaction. For simplicity, we assume the same 



trapping frequencies for the molecule and for the ion: 
LJi = ui m = to, which allows to decouple the center- 
of-mass (COM) and relative-motion degrees of freedom. 
We neglect the effects of the inelastic scattering that 
can occur in the collision of an ion and a ground-state 
(X E, v — 0, J = 0) molecule, which are not important 
when the particles remain well separated in the course of 
the cooling process. The relative motion can be described 
by the single-channel Hamiltonian 

H iel = ^ + ^ 2 (r-d) 2 + V(r), (9) 

where we assume V(r) — —C^/r 4 and we introduce some 
quantum-defect parameter ip, accounting for the effects 
of the short-range potential. A similar modeling turned 
out to be very successful in the case of ion-atom collisions 
[27l [2H ] , where the interaction at large distances is gov- 
erned by the same dispersion potential V(r) — —Ci/r A . 
The quantum-defect parameter ip can be interpreted as 
the short-range phase accumulated by the wave function 
due to the short-range molecular potential. This inter- 
pretation follows from the behavior of the solution of the 
Schrodinger equation with the Hamiltonian ([9]) . Perform- 
ing the partial wave expansion 9(r) = Ri(r)Yi m (9, (j>), 
we obtain the following short-range behavior of the radial 
wave function Ri(r): 

R t (r) ~ sin (if/r + w), r < y/R*/k,R*. (10) 



Here hk is the relative momentum, R* = ^j2C^/h 2 
is some characteristic length scale for the molecule- 
ion interaction, and ipi is interpreted as a short-range 
phase. The asymptotic solution (ITUl) fulfills the radial 
Schrodinger equation neglecting the energy fi 2 fc 2 /(2/i) 
and the centrifugal barrier H 2 1(1 + 1)/ (2 fir 2 ) terms, which 
requires r <C y/ R* Jk, R* . In this regime the phase ipi de- 
termined by the short-range interaction is independent of 
energy and angular momentum: ipi(E) = tp. This allows 
to enclose the effects of the short-range forces into a sin- 
gle parameter p, which together with Eq. (1101) implies a 
boundary condition on the wave function at r — > 0. The 
fact that ([T0|) fulfills the Schrodinger equation for k = 
and I — allows to express the s-wave scattering length 
a in terms of </>: a = —R* cot <fi [13] ■ The parametrization 
in terms of ip is valid as long as the distance Rq, at which 
V(r) starts to deviate from the asymptotic r -4 law (e.g 
due to the exchange interaction), is much smaller than 
R* , which should be well fulfilled in the case of ion-atom 
and ion-molecule interactions. 

In the framework of the quantum-defect model we can 
diagonalize the Hamiltonian @ with the boundary con- 
dition (fT0|) for small r, obtaining a set of eigenenergies 
E n (d) . By comparing with the predictions of the coupled- 
oscillator model we can estimate the range of distances 
where the latter is still valid. Fig. [5] shows the energy 
spectrum as a function of the trap separation d obtained 
from the numerical diagonalization of @ for some ex- 
ample value of ip, together with a few lowest energy lev- 
els calculated from the coupled-oscillator approximation. 
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FIG. 2. Energy spectrum of the relative motion for a molecule 
and an ion confined in spherically symmetric traps versus the 
distance d between traps, calculated for some e xampl e pa- 
rameters: ip — —tv/8, and R* — 3.48?, with I = \Jh/ /ioj. The 
exact numerical results (black lines) are compared with the 
prediction of the coupled-oscillator Hamiltonian (|17[) (thick 
red dashed lines). The vertical blue dashed line shows the 
minimal distance at which the Hamiltonian can be expanded 
quadratically around the equilibrium points, whereas the 
green dash-dotted line shows the distance at which the barrier 
between particles becomes equal to the relative energy of the 
particles. 

The figure depicts also the distance at which the barrier 
between the particles becomes equal to the relative ki- 
netic energy, and the minimal distance d m i n correspond- 
ing to w™ ax . We observe quite a good agreement between 
the oscillator model and the quantum-defect calculation; 
only close to d m - m docs the oscillator model start to de- 
viate from the accurate results. At such distances ne- 
glecting terms higher than quadratic in the expansion of 
the potential ceases to be valid. Apart from this, the 
two approaches differ in the vicinity of avoided crossings, 
representing trap-induced shape resonances [13, [H| be- 
tween molecular and trap states, which obviously cannot 
be incorporated into the coupled-oscillator model. 

D. Parametric resonance 

In general, the trapping frequencies for molecules 
trapped in the wells of an optical lattice are much smaller 
than the trapping frequency of an ion confined in an RF 
trap: u>i ^> ui m . In this case the exchange of energy be- 
tween the two particles is only possible in the presence 
of a time-dependent external potential, which can pro- 
vide the difference in energies between the phonons. In 
our scheme we propose to utilize the periodic modulation 
of the trapping frequency: io 2 (t) = ui 2 (l + 2A cos(cjfi)), 
where A is the amplitude of the modulation and ujf de- 
notes its frequency. This can be achieved by applying 
an additional potential, which is quadratic in space, and 
changes periodically in time. Periodic modulation of the 
trapping frequency gives rise to the phenomenon of para- 
metric resonance, allowing for exchange of energies be- 



tween the particles. In the subsequent analysis we con- 
sider the case when the modulation is applied to the 
ion trap, however, the exchange of phonons can be also 
stimulated by modulating the trapping frequency for the 
molecule, or the distance between the traps. 

In the case of an ion stored in an RF trap, a periodic 
modulation of the trapping frequency can be realized by 
adding to the electric field a component varying with 
the frequency ujf. The stability of a single ion in such a 
potential is analyzed in Appendix [Cj In the presence of 
modulation of the trapping frequency for the ion uif (t) = 
Lof(l + 2Acos(o;ft)), the total Hamiltonian H takes the 
form 

H = H Z +H rn , (11) 

A 2 
H m ^ -hu l cos(uj{t)[c a {a + a jt ) + Cf j {l3 + ^)] . (12) 

Here, c a and cp are the coefficients given in Appendix 
iBl resulting from the transformation diagonalizing the 
initial Hamiltonian. Now we transform to the interaction 
picture with respect to H z , and sufficiently close to the 
resonance (|w2 — ujf — uji \ <C lo\,uj2) we apply the rotating- 
wave approximation (RWA), leaving in the transformed 
Hamiltonian H m = cxTp(iH z t/h)H m exp(—iH z t/h) only 
slowly varying terms: 

H m — AhujiCaCpla^ exp(i(ojf + u>i — uj2)t) + h.c.]. (13) 

Transforming back from the interaction picture and mov- 
ing to a rotating frame with respect to the second os- 
cillator H = exp(—iH2t/h)H exp(iH2t/ti), with H2 = 
huifa^a, we obtain 

H RWA = ft(a; 1 +a; f )a t a+fiw2/3 t /3+^ c (a t /3+a^ t ). (14) 

This Hamiltonian describes a system of two coupled, sta- 
tionary oscillators with the coupling frequency fl c = 
AhujiCaCp, which in the weak-coupling regime (w c <C 
uji,ui2, — wi|) can be approximated by 

n c « A^fk?-^]- 1 . (15) 

Introducing the operators T\ — af3^ + T2 — i(a^ — 
a' /3) and T3 — a^a — j3' ft, fulfilling the same commuta- 
tion relations as the Pauli spin operators 

[Ti,Tj] = 2ie ijk Tk, (16) 

one immediately shows that (fT4")l is equivalent to the 
Hamiltonian of a two-level system driven by an exter- 
nal field. In the subspaces with a constant total number 
of phonons: N — a^a + ft^ ft — const, the Hamiltonian 
(fT4")l can be written as 

ffRWA = ^-T3 + MVZi, (17) 

where A = u>2 — — and we omit c-number terms 
affecting only the evolution of the total phase. Exactly at 
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resonance (A = 0), the expectation value of the difference 
of phonon numbers (T z ) evolves according to 

(T z )(t) = (T y )(0) ahx(0(t)) + (T z )(0) cos(0(t)), (A = 0) 

(18) 

where 9(t) = 2 J Q drf2 c (r). Hence, application of the ir- 
pulse 2 J °° cItQ, c (t) = n allows for exchange of the num- 
ber of phonons between the molecule and the ion, and, 
for instance, for realization of the SWAP gate based on 
qubits represented by motional states. For a modulation 
frequency u>{ that is off-resonance, and in the situation 
when initially all the phonon population is in one of the 
oscillators: (T z }(0) = N, (T x )(0) = (T y )(0) = 0, the 
time-dependence of (T z ) is governed by 



(T z )(t) = N 



1 



8tt 2 



tyj4£l 2 c + A 2 



4fi? + A 2 



(19) 

where we assume that fi c (*) = const. In this case the 
exchange of populations between the two oscillators is 
not complete, and one can observe that the decrease in 
the fidelity exchange has a typical Lorentzian shape with 
resonance width given by 2f2 c . 

So far we have discussed the results obtained in RWA, 
which are valid for A <C 1. In general, for larger A we 
can solve the time-dependent problem starting from the 
differential equations for the position operators in the 
Heisenberg picture derived from the Hamiltonian @ 



X + BX + 2 cos(w/i)CX = 0, 



(20) 



with the vector operator X T (£) = (xi(t) — Xi,x m (t) — x m ), 
and 



B = 



u; 2 2w c w 4 \ / Aujf , ( 



The form of Eq. (f2"U|) resembles a Mathieu differential 
equation, and analogously its solution can be found using 
Floquet expansion X = e lXt Y^=-oo v n e mulf *• Insert- 
ing this ansatz into (|2T)|) we find the recurrence relation 
for the operators v n : C(v n+ i + v n _i) = D n v n , where 
D n = I(\ + nuif) 2 — B, and I is the unitary matrix. 
The recurrence relation can be solved in terms of contin- 
ued fractions: v n = iJ n v„_i, v_„ = H- n v- n +i, where 
H n = [D n — CH n+ i}~ 1 C, H- ri = [D_ n — CH_ n _i)~ 1 C , 
and n > 0. The eigenfrequencies A have to be determined 
from the equation [C (.Hi (A) + -ff-i(A)) - D ] v = 0. 
For A <C 1, the lowest order approximation is obtained 
by assuming H±2 = 0, which leads to the formula (p~5|) 
for tt c and gives the following resonance condition: 



UJ2 — ^1 — Ldf 



0. 



(22) 



In comparison to the results obtained in RWA, for 
nonzero A, the resonance is shifted by the presence of 
the modulated trapping potential for the ion. This ef- 
fect is illustrated on Fig. 131(a), showing the position of 
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FIG. 3. Panel (a): modulation frequency uji that leads to the 
resonant coupling between the molecule and the ion (solid 
line with dots), and width of the resonance (dashed lines), 
calculated numerically from the continued fraction solution, 
compared with the analytical results (I22|l (solid line). Panel 
(b): effective coupling frequency Q c calculated numerically, 
compared with the analytical result of Eq. (|15|) (solid line). 
For both panels the calculations are performed with uj m = 
O.lwi and ui c = O.Olcji. 



the resonance as a function of A. It compares the ap- 
proximate result of Eq. (|2"2")) with the full solution of ([2H)) 
expressed in terms of continued fractions. In addition, 
Fig. [31(a) presents the width of the resonance Aoj, calcu- 
lated by assuming a Lorentzian shape for the resonance 
peak [3fj|. We observe that the resonance is broader for 
larger values of the amplitude A, however for A > 0.3 
its width does not further increase. The coupling fre- 
quency fl c as a function of A is presented in Fig. [31(b). 
This shows the analytical result (fT5)) and the exact result 
calculated from the continued fractions solution. We ob- 
serve that for small A, fl c depends linearly on A, and it 
is well described by the analytical formula. 



III. MOLECULE COOLING VIA PARAMETRIC 
RESONANCE 

To utilize the molecule-ion interaction to cool the 
molecules, one can either consider collective laser-cooling 
in the coupled molecule-ion system, or apply laser cooling 
to the ion separately, and then with the help of modula- 
tion of the trapping frequency, exchange the energies of 
molecule and ion. In this section we focus on the latter 
scenario. We distinguish the following steps: (I) the ion 
is laser cooled; (II) the particles are transferred to the dis- 
tance allowing for efficient exchange of phonons between 
them; (III) the modulation of the trapping frequency is 
switched on and the energies of molecule and ion are ex- 
changed by applying a 7r-pulse; (IV) the molecule and 
the ion are separated. Fig. [4] presents some typical dy- 
namics of the cooling process. It shows the molecule and 
ion energies at different cooling steps, while consecutive 
cooling steps are separated by vertical lines. 

We first focus on the problem of molecule and/or ion 
transfer to the distances where the sympathetic cooling 
can be efficiently performed. In order to avoid addi- 
tional heating due to the trap motion, we assume that 
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FIG. 4. Typical dynamics of molecule and ion energies for 
the cooling scheme where the ion is cooled separately from 
the molecule. The roman number corresponds to the different 
cooling phases described in the text. 



the molecule and the ion traps are moved adiabatically, 
which requires (1) z~k(t)/h <C u>k(t) for k — a,i, (2) 
\tl}j(t)\ < U3i(t)u} c {t),U32(t)uj c (t) for j = 1,2. The for- 
mer condition is necessary for the suppression of excita- 
tions during the movement of the traps, while the latter 
excludes the possibility of transitions due to the time 
dependence of 0Ji(t) and u>2{t), which change with the 
distance d{t). In general the suppression of excitations 
due to the motion of the ion's trap can be achieved by 
fulfilling the weaker conditions 31 [ 



t 



to 



dr [c a (T)fi(T) + d a (r)f m (T)} e**« 
dr ^(r)/,(r) + ^(r)/ m (T)]e*W 
with c a , cp, d a , dp given in Appendix |BI and 



< 1, 
« 1, 



fi(t) = u>iZi(t)/k + 2uj c (t)z m (t)/l m , 
fm(t) {t)/l m + 2u (t)zi(jt)/U 

ct>k(t)= /drcj fe (r), k = 1,2, 
Jo 



(23) 
(24) 



(25) 
(26) 

(27) 



where the limits to and t\ of the integral refer to the 
beginning and the end of the transfer process. Similarly, 
a weaker condition can be derived for the suppression of 
the probability of excitation due to the time dependence 
of u>i(t) and u>2{t): 



1 wi(r) +6j 2 (t) 
dr ; — e 



<1, fc=l,2 

:(0l(T)+0 2 (r)) 



4w c (t) 



/' 

■J to 



4w c (t) 



< 1, 



< 1, 



(28) 
(29) 
(30) 



In the phase (iii) when the energies of the molecule and 
of the ion are exchanged the pulse shape A(t) should be 
optimized in order to obtain the best efficiency of the 
transfer, and to avoid excitations due to switching on 
and off the time-dependent modulation potential. 



We performed numerical simulations of the cooling 
process starting from different energies of the molecule 
before the cooling process, and assuming that the ion is 
laser-cooled to the ground state: (rij(0)) = 0. For sim- 
plicity we have assumed that the transfer process is ideal, 
i.e. the conditions (|23]) - pg]l and ([28 ]) -(|30 ]l are fulfilled. In 
this way, after the transfer phase, the state of the parti- 
cles is described by a density matrix p = p a ® pp, where 
p a = J2 n Pn\ n )( n \ an( i P/3 = |0)(0| correspond, respec- 
tively, to the state of the first and the second harmonic 
oscillator of the Hamiltonian ((7} . Fig. [5] shows the final 
mean phonon number of the molecule for u> m = O.lwi, 
uj c = O.Olwi. Fig. 03(a) shows the length t{ of the ex- 
change pulse as a function of the amplitude A, while 
the dashed line presents the lower bound for the length 
due to the physical coupling uj c between the oscillators: 
imin = t/(2w c )- F° r simplicity, in the calculations we 
have assumed a rectangular shape of the pulse for the 
amplitude of the parametric modulation A(t). Hence the 
pulse length is given by t- lmp — 7r/(2f2 c ). The sudden 
turning on and off of the parametric modulation leads to 
an additional heating that increases with the amplitude 
A. This is shown in Fig. 0(b), presenting the final mean 
phonon number of the molecule (n m ) for different values 
of the initial mean phonon number (n m (0)). We observe 
that, even without using an optimized pulse for A(t), the 
heating introduced by imperfections of the transfer pro- 
cess remains quite small for A < 0.3, and is practically in- 
dependent of (n m (0)). We notice that the increase of the 
molecule energy with A is not monotonic and that there 
is a minimum in the molecule energy around A w 0.2. On 
the contrary, the molecule and ion energies for detunings 
far from the parametric resonance, when the exchange of 
energy is suppressed, increase monotonically with A. 



IV. SIMULTANEOUS LASER COOLING AND 
ENERGY EXCHANGE 

We turn now to the description of laser cooling in a 
coupled system consisting of a trapped molecule and an 
ion, where laser cooling acts at the same time as the para- 
metric modulation. This process consists of three sepa- 
rate phases, (I) transfer of particles to distances allowing 
for efficient exchange of phonons between the molecule 
and the ion; (II) laser cooling and simultaneous exchange 
of energy through parametric resonance; (III) separation 
of the particles. 



A. Master Equation 

We assume that the cooling laser affects only the ion, 
and that the coupling between molecule and ion is suffi- 
ciently weak: u> c -C 0,7, where Q is the Rabi frequency 
due to the interaction of the ion with the laser field, and 7 
is the spontaneous emission rate. In this way, the absorp- 
tion and spontaneous emission processes of the ion are 
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FIG. 5. Panel (a): Length of the rectangular 7r-pulse (solid 
line), and lower bound for the pulse length (dashed line) due 
to the physical coupling ui c . Panel (b): Final mean phonon 
number of the molecule (n m ) as a function of the amplitude 
A and of the initial mean phonon number of the molecule 
(n m (0)). Calculations are performed with ui m = O.ltJi, ui c = 
O.OlWi, and (n*(0)) = 0. 



not affected by the coupling with the molecule. Similarly 
to the considerations presented in the previous sections, 
we restrict ourselves to the case of one-dimensional mo- 
tion along the direction determined by the centers of the 
traps. Using the standard theory of quantum stochastic 
processes [H, [H| , we derive a master equation (ME) for 
the reduced density operator p obtained by tracing over 
the empty modes of the radiation field. In the frame 
rotating at the laser frequency u>l, the ME reads J34| 

9 =\ [H + H lnt + H las ,p] 
in 



duW{u)e lkz ' u pe- lkz ' u a+ (31) 



(7 + (T_p — p<7+(7- 



Here, H denotes the Hamiltonian (fTTj) for the external 
degrees of freedom, Hi nt is the Hamiltonian for the in- 
ternal degrees of freedom 



Hi, 



hA 



(32) 



where A = ujl—luq denotes detuning, ujq is the frequency 
of the optical transition, and a z , a± are the standard op- 
erators describing transitions in the two-level molecule. 
The Hamiltonian Hi as describes the interaction of the ion 
with the laser field 



Ml 
~2~ 



e '<t + 



h.c. 



(33) 



with k denoting the wave vector of the laser field. Finally 
W(it) is the angular pattern of the spontaneous emission, 
which for dipolar transitions is W = 1(1 + u 2 ). 

We focus now on cooling in the Lamb-Dicke regime, 
where the laser wavelength is much larger than the ion's 
trap size. The presence of the small Lamb-Dicke pa- 
rameter r\ = k\Jhj '(2miU>i), allows to expand the ME 
in powers of r\. To lowest order in r\ the evolution of 
the internal and external degrees of freedom is decou- 
pled, and the internal state can be adiabatically elimi- 
nated (see Appendix [Dj . The adiabatic elimination is 
valid provided that the internal dynamics is faster than 
the lowest-order coupling between the internal and ex- 
ternal degrees of freedom, which requires rfil -C 7, Q. As 
a result of the adiabatic elimination, we obtain the fol- 
lowing ME for the reduced density operator p e obtained 
by tracing over the internal degrees of freedom 



9e =~r 

in 



1 T Ml ,,22 

H - — KVylsKz-i - —{a x ) s k Zi ,p e 

+ ^rk 2 (a + a-) s (2z l p e z i - zf p e - p e zf) 



2n2 r°° 



k 2 n 



(zi(T)p e Zi - ZiZi{r)p e ) 
x ((ff,(T)ff v (0)}, - (a 



2 

ills 



h.c 

(34) 



where we include terms up to the second order in 77. Here, 



Zi(r) = e TCoE Zi = e- UH / h Zl e ltH/r \ 
(<r y (T)(jy(Q)) s = Tr I {a y e TCol a y p ss } 



(35) 
(36) 



and (,) s = Tr i{.p s } denotes an average taken in the sta- 
tionary state of the density matrix p s of the internal de- 
grees of freedom. It is defined by £oi9s — 0, where the 
Liouvillians Cqi and Cqe are given by 



£019 : 



1 



-(a- +<7+) 



hA 



2<j-pa+ — a+cT-p — po+<j-\ , 



£oeP =tt [H, p] 

ih 



(37) 
(38) 



The two-time correlation functions (a y (T)a y (0)} ss can be 
calculated with the help of the quantum regression the- 
orem (see Appendix [E] for details) . The first term of the 
ME (|34|) describes the evolution due to the trapping po- 
tentials and the coupling between the molecule and the 
ion. We observe that the trapping potential for the ion is 
modified by the interaction with the laser field. The sec- 
ond term of (|34p represents damping due to spontaneous 
emission, while the third term describes the excitation of 
the ion due the cooling laser. 
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B. Results in RWA 



where 



In this section we discuss laser cooling in the molecule- 
ion system on the level of RWA. Although a treatment 
based on RWA is in principle limited to the regime A <C 
1, the analysis presented in this section provides a good 
qualitative description of the cooling dynamics, and it 
can be carried out fully analytically. A description of the 
laser cooling for arbitrary A requires going beyond RWA, 
and this is addressed in the next section. 

In the framework of RWA we make the following as- 
sumptions: (i) H is approximated by Hrwa given by 
(TT4"]) ; (ii) RWA is applied to the second and to the third 
term in the ME (|34[) . which results in omitting the 
quickly rotating terms; (iii) we neglect the second and 
the third term in the commutator in Q34p , describing the 
effects of the laser field on the ion's trap, which requires 
?yil Cwj. In this way the ME (l34l) reduces to 

• _ 1 r ~ 

in L 

+ {s-{ap e fi - fiap e ) + S+(a^p e (3 - pet p e ) 
+ (S+ + D){fip e p - Pe ) + /i.e.}, 

(39) 

where D — 7] 2 j(<j+<j^) s /5, 

v 2 n 2 r°° 



S(w) = 



S v 

s,. 



4 Jo 

v 2 n 2 r°° 



drf*(r) ((a y (T)a y (0)) s - (a y ) 2 ) , (40) 

> 

drf v (r)((a y (r)a y (0)} s -(a v } 2 a ), (41) 



for v = a,b. Here f v (r) are the functions determined by 
the evolution of the operator j3 in the Heisenberg picture 

P(t) = e it&RWA/n 0e- itSRWA/h = f m (r)a+f b (r)0. (42) 

We observe that the ME (f3T)| contains, apart from the 
terms describing the laser cooling of ion, also mixed terms 
involving a and ft operators corresponding to the action 
of the laser field on the molecule. In the considered 
regime of weak coupling w c <C 7, fi, mixed terms can 
be neglected in comparison to the terms describing the 
laser cooling of the ion only. This can be verified by in- 
serting the solutions for / m (r) and /&(r), which exactly 
at parametric resonance + lo{ =1-02) takes the form 
/ m (r) = ~ie-^ 2t sin(fi c i) and f b (r) = e^" 2 * cos(Q c t). 
Inserting these solutions into (j40 |) - (|4T]) we obtain 



S~ 



I [Sfa + do) - 5( W2 - n c )] , 
\ [s{lu 2 + n c ) + s(uj 2 - n c )] , 
\ [s(-u 2 - n c ) - s(-u2 + n c )} , 
\ [S(-ui - n c ) + s(-oj 2 + n c )} , 



(43) 
(44) 
(45) 
(46) 



dre™ (K(rK(0)) s -K)2), (47) 



The function S(uj) is calculated in Appendix [El and 
from the form of this solution we observe that S(to) is 
determined by the laser cooling parameters 7, f2, and 
A. Therefore in the regime f2 c -C 7, f2, A we have 
S(u)2 + f2 c ) ps S(ui2 — O c ) ps 5(^2) and 5(— W2 + ^ c ) ps 
S(— u>2 + dc) ~ S(— W2). Hence, in the lowest approxi- 
mation the coefficients 5^" and becomes independent 
of the coupling between molecule and ion, and identi- 
cal with the coefficients describing the laser cooling of 
a single ion. On the contrary, in this regime the coeffi- 
cients S~ and S£ become very small, therefore the terms 
describing the damping process and involving a and 
operators can be neglected. 

From the ME (f3"9"]l we derive a set of differential equa- 
tions describing the dynamics of the cooling process. 
This straightforward calculations leads to 

it = iujv + (a + — a-)u, (48) 

v = iuju + (a_|_ — a,— )v — i2uj c (w — s), (49) 

s = iuJ c v + 2(a + — a_)s, (50) 

w = —iui c v, (51) 

where the variables u, v, s and w are defined by s = 
08t/J) - (p0)„ w = (at a) - (/?t/3) s , u = C8at) + (^ a ), 
v = (/3t a ) - (/3 a t). Here (/3t/3) fl = a+/(a_ - a+) de- 
notes the final energy of the ion at the end of the cool- 
ing process, a± = Re(S* + D), lo± = lm(S b + D) and 
lu = u>2 — — + uj + + w_ . We observe that the final 
(equilibrium) energy of the molecule is the same as for 
the ion in the absence of the molecule (34j. The set of 
differential equations (l48l) - (|5Tj) has one pair of complex 
eigenvalues and one pair of real eigenvalues. One can 
easily verify that the larger of the two real eigenvalues 
provides the cooling rate of the molecule 



r 

2 



1 

71 



■4toi 



in 



uj 2 T 2 



1/2 



Ir 2 



(52) 



where T = 2(a_ — a+) is the cooling rate of the ion in 
the absence of the molecule. At the parametric reso- 
nance we have to = u + + cj_, and for typical parame- 
ters corresponding to the sideband cooling regime we get 
cj + ,u;_ <C r. The cooling rate of the molecule T m as 
a function of the cooling rate of the ion T, for different 
values of the parameter u> is shown in Fig. [5] We note 
that for r < 4J7 C the molecule is cooled with a rate T/2. 
Further increasing of T does not result in faster cooling of 
the molecule, and in the regime V > 4fi c the cooling pro- 
cess becomes ineffective. This behavior is also illustrated 
in Fig. showing the cooling rate versus the amplitude 
A of the parametric modulation, calculated for different 
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FIG. 6. Cooling rate of the molecule F m as a function of the 
cooling rate of the ion F, for different values of the parame- 
ter uj (see text for details). Cooling rates are scaled by the 
effective coupling Q c . 
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FIG. 8. Dynamics of the mean phonon number of the 
molecule n mo i(t) and the ion rnoaif)- The calculations are 
performed for the parameters of Fig. [7] and Q. = Q.luji. The 
dotted horizontal line shows the final energy for the molecule 
and the ion. 
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FIG. 7. Cooling rate of the molecule Y m versus the am- 
plitude of the parametric modulation A, for different values 
of the Rabi frequency Q. For comparison we include a line 
showing the values of 2Q C (solid line without symbols), and 
the cooling rate of the ion r (horizontal lines). The results 
are obtained for cj m = 0.2oJi and uj c = O.Oluji, m m = m; and 
for cooling parameters corresponding to the sideband cooling 
regime: 7 = O.lcji, A = — 0Ji and r\ — 0.1. 



values of the Rabi frequency fl. The presented results 
are obtained for the following values of the parameters: 
LJ m = 0.2^, lo c = O.Olwj, m m = im, 7 = O.la;,, A = -w, 
and rj = 0.1. For comparison we include a line showing 
the values of 2f2 c . We observe that for too small values 
of A (when 4£l c < T) the cooling is inefficient and the 
cooling rate of the molecule is smaller both than F and 
£l c . For larger values of A, when Afl c > T, the cooling 
rate T m becomes independent of A and equal to T/2. 



C. Beyond RWA 



In this section we consider cooling in the ion-molecule 
system taking into account effects beyond RWA. We start 
from the ME (|34[) with H represented in terms of momen- 



tum and position operators 
1 



2 

Pi/ 1 1 2 2 



2huj c 



(53) 



AmiUJ i z i cos(wfi), 



and we derive equations of motion for average values 
quadratic in the operators z t , pi, z m , p m . In principle 
the term linear in z% in the commutator in (1341) cou- 
ples the observables quadratic in positions or momenta 
to the linear ones. However, similarly to the analysis per- 
formed in RWA, we can neglect it assuming the effects of 
the laser field on the trapping potential for the ion are 
small for rjU, <C Wj. In this way we obtain a set of lin- 
ear equations, which we discuss in Appendix |FJ We solve 
the equations of motion numerically, determining the dy- 
namics of the cooling process and the final energies of 
the molecule and the ion. Some sample dynamics of the 
cooling process is presented in Fig. HI which depicts the 
mean phonon number of the molecule and of the ion as 
a function of time. The presented results are calculated 
for uj m = 0.2cj;, cu c = O.Olwi, m m = wii, and for cool- 
ing parameters corresponding to the regime of sideband 
cooling: 7 = = O.lw,, A = — uj; L and 77 = 0.1. Fig.[5](a) 
presents the dynamics for A = 0.01, for which T > 4fi c 
(cf . Fig. [7]) . In this case we observe a very weak exchange 
of energies between the particles, which is due to the too 
fast cooling of the ion. This behavior is somehow similar 
to the Zeno effect that can be observed in quantum sys- 
tems, corresponding to the suppression of the dynamics 
in the presence of frequent measurements performed on 
the system. In our case the dynamics of energy exchange 
is suppressed due to the too strong damping of the ion's 
motion. On the other hand, for A = 0.1 the parameters 
are well within the regime of efficient cooling (T < 4Q C ), 
which can be deduced from Fig. [7] In this case the dy- 
namics exhibits the presence of oscillations, which are 
damped in time due to the laser cooling process, which 
is illustrated in Fig. [SI (b). 

Fig. [5] shows the dependence of the final energy of the 
molecule and the ion on the amplitude A, calculated for 
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FIG. 9. The final mean phonon number of the molecule and 
the ion as a function of the amplitude A. The calculations 
are performed for the same parameters as in Fig. [8] 




FIG. 10. The final energy of the ion as a function of the 
detuning A, and for different values of the amplitude A. The 
results are calculated for the parameters of Fig. [Jj 



the same parameters as in Fig. [3] We observe that the 
final energies of the molecule and the ion are equal, and 
that they increase with the amplitude A. Hence, the pres- 
ence of the oscillating potential decreases the efficiency of 
laser cooling. To investigate whether this effect is related 
with the coupling between the molecule and the ion, or 
if it originates purely from the parametric modulation 
of the ion's trap, we performed simulations of the cool- 
ing, putting w c = 0, which describes the situation of the 
molecule separated from the ion. 

Fig. [TU] shows the dependence of the final energy of 
the ion on the detuning A, calculated for different val- 
ues of the amplitude A. We observe that the modulation 
of the trapping frequency leads to an increase in the fi- 
nal energy, similarly to what is observed in the coupled 
ion-molecule system. For the parameters used in the sim- 
ulation, the cooling is performed in the resolved-sideband 
regime, and for non-zero values of A we observe the ap- 
pearance of a second sideband at the frequency u>i + ujf. 
For higher values of A (A ~ 0.1), the cooling tuned to 
this second sideband is almost as efficient as the cooling 
at the sideband located at cjj. 

Figs. QT] and [12] shows the final energy of the ion as a 
function of the Rabi frequency and of the spontaneous 
emission rate, respectively, evaluated for different values 
of the amplitude A. Again, we note that for any combi- 
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FIG. 11. The final energy of the ion as a function of the Rabi 
frequency Q, and for different values of the amplitude A. The 
results are calculated for the parameters of Fig. [JJ 




FIG. 12. The final energy of ion as a function of the spon- 
taneous emission rate 7, and for different values of the am- 
plitude A. The results are calculated for the parameters of 

Fig.m 



nation of the cooling parameters the final energy of the 
ion increases with A. Finally, we can conclude that the 
lowest energies are achieved when f2, 7 <C u>%, and for the 
laser tuned to the first red-detuned sideband A = — w,-. 
At the same time, the amplitude A should be kept as 
small as possible, limited only by desired rate of the en- 
ergy transfer between the molecule and the ion. 



V. CONCLUSIONS 

In summary, we have analyzed sympathetic cooling of 
the center-of-mass motion of a polar molecule prepared 
in the internal rovibrational ground state and interact- 
ing with a laser-cooled ion. The ion and the molecule are 
confined in separate trapping potentials, and the energy 
exchange is mediated via the long-range ion-molecule po- 
larization potential and the additional ion trap modula- 
tion allowing for an efficient exchange of phonons be- 
tween the two particles. We have considered two dis- 
tinct cooling schemes. In the first method the laser 
cooling and the energy exchange pulses are separated, 
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while in the second scheme these two processes happen 
at the same time. In both cases the resulting cool- 
ing rate is determined by the slower of these two pro- 
cesses, which typically is the energy exchange process 
stimulated by the parametric resonance effect. For an 
example system of a Ca + ion interacting with a KRb 
polar molecule, the cooling rate of the ion can be as 
large as I\ = lOOKHz 35], while the maximal coupling 
frequency uj™ ax determined in Section [TT] is limited to 
2ir x 27KHz for trapping frequencies u>i — 2tt x 1MHz 
and oj m — 2tt x 100kHz. Assuming a slightly smaller 
value for the coupling ui c = 2tt x IOKHz, for which the 
linear expansion of the potential is applicable, one gets 
the maximal cooling rate T = 2fi c ps 6kHz for a mod- 
erate value of the modulation amplitude A — 0.3. This 
can be further increased up to V ss 2uj c = 20kHz for a 
modulation amplitude Apz 1 , at the expense of a higher 
energy of the molecule at the end of the cooling process. 
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Appendix A: Derivation of the effective Hamiltonian 
for trapped ion and molecule 

In this Appendix we present a microscopic derivation 
of the effective Hamiltonian (QJ in an adiabatic approx- 
imation and we discuss its validity. We consider an ion 
of mass m% with positive charge +e. For simplicity we 
do not take into account the internal structure of the 
ion, which has a different electronic structure than the 
molecule, and is not affected by the optical trapping po- 
tential. We assume that the diatomic molecule is in its 
ground vibrational state (v = 0), and we treat it as a 
rigid rotor, neglecting the couplings to excited vibrational 
states. The molecule consists of two nuclei denoted as A 
and B, with charges Z A e and Z^e, and masses ttla and 
77i b , respectively. For simplicity we carry out our deriva- 
tion in the nonrelativistic approximation, neglecting the 
effects of fine and hyperfine level splitting. The total 
Hamiltonian can be written as 



H = Hi 



kin 



H e \ + H\ as + H T f + Hi, 



The kinetic energy of the nuclei is 

.1? r.2 



Pi 

2m,- 



Pr 

2M 



BJ 2 



(Al) 



(A2) 



where the labels i and m denote, respectively, the ion 
and molecule COM degrees of freedom, M = ra A + t>ib, 
BJ 2 is the kinetic energy of a rigid rotor, while B de- 
notes the rotational constant of the molecule in the vi- 
brational ground state (v = 0). In our derivation all 



the unprimed quantities refer to the space-fixed frame, 
whereas the primed ones refer to the rotating, molecule- 
fixed frame. 

The term H e \ is the electronic Hamiltonian, including 
the kinetic energy and Coulomb interactions 



E 



p±_ 

2m r 



Z A c 2 



3>j' 



(A3) 

with the label j enumerating the electrons of the 
molecule. 

The term _ffi as describes the interaction of the molecule 
with a far-off-resonance laser light creating an optical 
potential. In the electric dipole representation 



Hus = — dEi as (r m , t), 



(A4) 



where d = e(Z A r A + ZBVB — Ylj r j) ^ s tne dipole moment 
of the molecule. We have applied the long-wavelength 
approximation, neglecting changes of the electric field 
Ei as on the scale of the molecule. The molecule can 
be optically trapped, for instance, by a pair of counter- 
propagating laser beams with circular polarization [36| . 
creating a standing wave 

E las (r,7j) = cos(fc L a; + <j> x ) 2 (E e+e lULt + c.c.) 

+ sm(k' L x + (j) x ) 2 (E'oe+e^ + c.c.) 

+ (x y z). (A5) 

Here, the abbreviation (x — > y — > z) denotes the rest of 
the terms obtained by cyclic permutations of the x, y and 
z coordinates. For each Cartesian direction {e x ,e y ,e z } 
we include two pairs of laser beams with amplitudes Eq, 
E'q, wave vectors k^, k' L , and frequencies ujl = k^c, 
lu' l = k' L c. This allows to cancel the influence of the 
tensor shifts, inducing a position dependent splitting of 
the rotational levels [36| . For simplicity we have assumed 
the same amplitude, wave vectors, and frequencies for 
all three directions of the laser beams. The versors of 
the spherical basis are denoted by {e_,eo,e + }, where 
eo = e z , e± = T( e x ± i^y)/V^i and <t> x ,y,z are the phase 
factors. 

The Hamiltonian of the radio-frequency (RF) field if r f 
creating the trapping potential for the ion reads 



H ri = e$ rf (r l ,t) - dE rf (r m ,i), 



(A6) 



where $ r f(r,t) is the time-dependent electric field of the 
RF trap: 



1 



5> r f(r,t) =—(u x x + u y y +u z z )cosw r ft. 



(A7) 



Here, w r f is the frequency of the time-dependent electric 
potential, and Uk (k = x, y. z) are amplitudes depend- 
ing on the trap geometry [371 ] . The electric field at every 
instant of time fulfills the Laplace equation A$ = 0. 
Hence, the coefficients u^, Vk are subject to the follow- 
ing condition: u x + u y + u z = 0. For simplicity we have 
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assumed that the RF trap is created only by the time de- 
pendent component, which allows to avoid the effects of 
molecule polarization in the electric field of the RF trap. 
In the second term of (|A6I) describing the interaction of 
the molecule with the RF field, we have assumed the size 
of the molecule to be much smaller than the characteristic 
length scale of the electric field E r f(r, t) — — V$ r f (r, t). 

The last term in the Hamiltonian describes the inter- 
action of the molecule with the ion 



-dEi(r r , 



(A8) 



where E,(r) = er/r 3 . We have only included the lowest 
order term in the multipole expansion, assuming the ion- 
molecule distance to be much larger than the size of the 
molecule. This breaks down at shorter distances, when 
the ion starts to distinguish separate components of the 
molecule, however, the details of the short-range physics 
of the ion-molecule interactions are not relevant for our 
analysis. 

Below we indicate the basic steps of the derivation. 

Expansion in the basis of Born-Oppenheimer wave 
functions for the electron motion. We start from gen- 
erating a complete set of electronic wave functions, pa- 

I 



rameterized by the positions of the molecule COM, r„ 
and by the orientation (9, </>) 



-ff i$n({rj}|r n 



= E n $ n ({T j }\r m ,0,cf>), (A9) 



where {r^} denotes the set of electronic coordinates. In 
this way the total wave function can be expanded in the 
basis of Born-Oppenheimer electronic wave functions 



n 

$„(fo}|r m ,M). 



(A10) 



Since the basis is complete, the expansion of the wave 
function does not involve any approximations. 

Born-Oppenheimer approximation. We use the Born- 
Oppenheimer method, treating the electron motion in the 
adiabatic approximation. This assumes that the time 
scale of the electron dynamics is much faster than the 
dynamics of the molecule's nuclei, which is typically ful- 
filled since the electron is much lighter than the other two 
particles. Substituting expansion (|A10|) into the time- 
dependent Schrodinger equation, we obtain a set of cou- 
pled differential equations for the expansion coefficients 
c n (ri,r m ,6,(f>,t): 



ih 



dCn 

dt 



ffkin + e^rfOi, t) + E n 



$„|dE 

las ipm 1 

t) + dE r f(r m ,t) + dEj(r m - rj)|$ r . 



(All) 



r 



Adiabatic elimination of the excited electronic states: 
derivation of the optical trap potential. For a far-detuned 
laser, the excited states are only weakly populated and 
can be adiabatically eliminated. In the spirit of the 
Born-Oppenheimer approximation we can assume that 
the transitions between ground and excited electronic 
states, due to the laser light, occur on a time scale much 
shorter than the motion of the molecule and the ion, and 
the dynamics of their COM motion can be decoupled 
from the internal dynamics. The result of the adiabatic 
elimination is most conveniently expressed in terms of 
the dynamic polarizabilities ai|(a>) and aj_(io) in the di- 
rections parallel and perpendicular to the internuclear 
axis [38[. The calculation can be performed in a similar 
manner as described in Ref . [36j , with the only difference 
that here we consider a three-dimensional optical lattice. 
As a result we obtain the following trapping potential for 
a molecule: 



Kpt(r) = -V a + \Mui 2 m {v-ry 



(A12) 



where Vq is the position independent component of the 
AC Stark shift V = \E \ 2 [§a|| (w L ) + %a±(w L )] , and 



the trapping frequency is given by 
2 



M 



\E k L \ 2 [fa||(wi) + 



\E' k' L \' fa,|K) + 



(A13) 



The second laser with frequency ui' L is blue detuned from 
the resonance in order to make and aj_(w) both 

negative. With this assumption uj 2 n is positive, and one 
can additionally adjust the parameters of the second laser 
to fulfill 



\E k L 



a \\ 



K^| 2 [fa„K)-|a ± K)] 



0, (A14) 



which allows to cancel position-dependent terms sensitive 
to the orientation of the molecule. We note that in the 
case of a three-dimensional lattice created by laser beams 
of identical intensity, the AC Stark shift in the center of 
the lattice well (r = 0) does not depend on the molecule 
orientation, in contrast to a one-dimensional optical lat- 
tice m. 

Time averaging over fast oscillations of the RF field: 
derivation of the Paul trapping potential. We replace the 
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time-dependent RF field by an effective adiabatic poten- 
tial for the ion, neglecting its fast micromotion on the 
time scale of tt> r f [37j |. Following the standard approach 
described in [33], we obtain 



Vrf (r) = -m l (u x x 2 + u y y 2 + lu z z 2 ) 



(A15) 



where Wk = (q 2 ./2) 1 / 2 u I f /2, with qu = 2euk / (m2W r 2 f ) for 
k = x,y,z. We can estimate the influence of the RF 
field on the molecule by considering the corresponding 
AC Stark shift fi 2 /A, resulting from the transitions J = 
— » J = 1. For a typical molecule (d ~ ID) and for 
electric fields down to distances of the order of few tens of 
harmonic oscillator lengths from the trap center, the Rabi 
frequency ft <C w r f. The rotational splitting given by B 
is typically much larger than uj r f < 10MHz, therefore 
the detuning is of the order of A ~ B ~ 1GHz and the 
resulting AC Stark shift of the molecule in the RF field 
can be neglected in comparison to the frequency of the 
optical potential. 

Now, including the adiabatic elimination and time- 
averaging of the RF potential, we obtain the following 
equation for the wave function c g (i"i, r m , i) corre- 
sponding to the electronic ground state: 



dt 



H kin + V r{ (Ti) + V opt {r m ) + E g 



dgg (r m r i ) 



Cg(Y l 



?,0,t). (A16) 



Here, d gg = ($ 5 |d|<i> g ) is the dipole moment of the 
molecule calculated in the electronic ground state. It 
depends on the orientation of the internuclear axis, and 
its coordinates in the space- fixed spherical basis are given 
by ( d gg)q = d gg ■ q = d gg Cf' (6,<j>) for q = 0,±1. 

Adiabatic elimination of the molecule rotational degrees 
of freedom. In analogy to the Born-Oppenheimer approx- 
imation for the electronic wave functions, we can perform 
an adiabatic approximation with respect to the ionic and 
molecular translational degrees of freedom. This can be 
done since the molecular rotations occur on a timescale 
2ir/B, which is much shorter that the time scale of the 
translational motion, given by 2ti juii and 2n /u> m for the 
ion and the molecule, respectively. The total wave func- 
tion can be expanded in the basis of wave functions 
A„(#, 4>\vi, r m ) of the rotational motion, parameterized 
by the positions of the ion and of the molecule 

c g {Ti,T m ,6,4>,t) = ip n {rj, r m , t)A n (9, <f)\rj, r m ), 



which fulfill 



(A17) 



[BJ 2 - d flff E i (r m - ri)] A„((9, ^r^r* 



= £; n (r m -r i )A„(e,^|r i ,r m ). (A18) 
The solutions of (|A18[) can be found by expanding in the 



basis of spherical harmonics Yjm(&, 4>) 



A,, 



jmYjm 



(A19) 



JM 



where the z axis is chosen along r = r m — r^. Introduc- 
ing a dimensionless parameter r\ = d gg Ei(\r\)/ B charac- 
terizing the coupling between angular momenta J and 
J ± 1, we observe that typically r\ <C 1 at ion- molecule 
distances corresponding to our cooling scheme. There- 
fore the eigenstates A„ are dominated by | J, M) with a 
small admixture of |J± 1,M) proportional to rj. In this 
case we can solve for E n and A„, taking into account 
only | J, M), |J± 1,M) states, and neglecting couplings 
to other states. This yields 



Ejm — B 



J(J+1) + 



J(J+ 1) - 3M 2 



2 J{J + 1)(2J- l)(2J + 3) 

(A20) 

which is the sum of the unperturbed eigenvalue plus a 
small correction. For a molecule in the rotational ground 
state the eigenenergy is E 00 = —Brj 2 /6, which leads to 
an attractive potential between the ion and the molecule 



V int (r) 



(dggef 1 

6B r 4 



(A21) 



Appendix B: diagonalization of the coupled 
oscillator Hamiltonian 



The Hamiltonian ([5]) can be diagonalized with the help 
of the following transformation: 



a f + a = c a (a 1 + a) + c (ft + (3), 
b> + b = d a (a^ +a) + dp{fi + /?), 



6+ -b = W[c a (ft ■ 



a 



0) - c ( g(a t - a)] , 



(Bl) 
(B2) 
(B3) 
(B4) 



where 



luJi(u} 2 - LU 2 ) 



cp 



Iwiiw 2 - ui 2 ) 



, (B5) 



luj m (u} 2 - Ujf ) 



luj m (u} 2 - to 2 ) 



n 2 



LU 2 



, (B6) 



u>i and Cl>2 are the frequencies of the two eigenmodes 
cj 2 2 = \ [uj 2 + u)^ =F with the upper (lower) sign re- 
ferring to the frequency u)\ (w 2 ), respectively, and 

& = ^ 2 - c^J 2 + 16 W 2 w m ^, (B7) 

(B8) 



W = 
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Appendix C: periodic modulation of the trapping 
frequency for an ion confined in an RF trap 



For an ion stored in RF trap, a periodic modulation 
of the trapping frequency can be realized by adding to 
the electric field an additional component, which varies 
with a frequency uji. To be more precise, the electric 
potential felt by the ion in the center of the trap should 
be of the form $(r) = (U + V cos uj r ft+W cos uift)(a x x 2 + 
a y y 2 + a z z 2 ) / (2tq), where w r f is the frequency of the RF 
field, ro is the characteristic distance describing the size 
of the RF trap, and a x , a v , a z are geometrical factors 
depending on the spatial configuration of the trap [37| . 
At every instant in time $(r) has to fulfill the Laplace 
equation: A<£> = 0, hence a x + a y + a z = 0. Classically, 
the motion of the ion in the direction k (k = x, y, z) is 
governed by 



d 2 x k 



'o 



■Vcosuirit + Wcosui{t)x = (CI) 



For W = 0, Eq. (|C1|) reduces to the well known Math- 
ieu equation. Its lowest stability region is presented in 
the first panel of Fig. [13] as a function of the parame- 
ters a — AaueU / {miW 2 ^) and q — 2aueV '/ '{rniU) 2 { r 2 ) . 
The motion is stable in the gray-shaded regions of the 
a-q plane. Inclusion of a field varying with frequency u>{ 
modifies this stability diagram, which can be seen from 
Fig. [13] presenting the regions of stability calculated for 
Wf/u; r f = 0.05 and A = 0.1, 0.2, and 0.3. At each point 
on the a-q plane, the amplitude W is adjusted in such 
a way that the motion of the ion, averaged over the fast 
micromotion on the time-scale of w r f, can be described 
effectively as a motion in the harmonic trap with the 
modulated frequency uj 2 (t) — cj? (1 + 2Acos(uift)). Fi- 
nally, the parameters of the RF trap have to be chosen 
in such a way that the motion of the ion is stable simul- 
taneously for all three spatial directions. 



Appendix D: Adiabatic elimination 



In this section we perform adiabatic elimination of the 
internal dynamics. The basic steps of this derivation 
are similar to the treatment of laser cooling of a sin- 
gle ion (cf. for instance Ref. [Hj]). First we expand 
the ME (f3"Tj) in powers of the Lamb-Dicke parameter 
rj = kyJh/(2m,iijJi) <C 1, retaining terms up to the second 
order in 77: 




FIG. 13. Stability diagram of Eq. ([Cl]) for w f /w rf = 0.05, 
and ,4 = 0, 0.1, 0.2 and 0.3. The amplitude W of the com- 
ponent varying with the frequency uj{ is determined from the 
parameters a, q and A (see text for details). 



where 
C01P - 



1 

ih 



(0) 



T ro 

2 Vv-pv 



Hint , P 

(T+CT-P - P<T+<T-} , 



£oeP = -r [H, p] 
in 



dp - 
C2P ■ 



1 

ih 
1 

ih 



n \ as > p 

n \ as > P 



(D2) 
(D3) 
(D4) 



^fc 2 CT_ {2z l pZ i - Z 2 p - pzf) <7 + . 



(D5) 



Here Hj® , H,^! , and are the terms of expansion of 
H\ as in the powers of 77: 



r(2) 



H 



H 



H 



(0) 

las 

(i) 

las 

(2) 
las 



(cr_ + C7+), 

(ikzi<j+ + h.c), 



- — ( fc 



h.c). 



(D6) 
(D7) 
(D8) 



p = (Cqi + C 0E + £1 + £2) p(t), 



(Dl) 



To lowest order in rj, the dynamics is described by terms 
£01 and £qe, and the evolution of internal and external 
degrees of freedom is decoupled. The adiabatic elimi- 
nation can be done provided that the internal dynamics 
is faster than the lowest-order coupling between the in- 
ternal and external degrees of freedom, which is given 
by £\. Since the time-scale of the internal dynamics is 
given by and 7 (cf. Appendix [Ej, the condition for the 
adiabatic elimination can be written as r/Q <C 7, £1. 
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The adiabatic elimination is done with the help of the 
projection operators technique [32| . We introduce the 
operator V p = p s Tii{p}, where the trace is taken over 
the internal degrees of freedom, and p s denotes the sta- 
tionary density matrix of the internal states: £aip s = 0. 
First we project the ME (|Dlj) using V and the comple- 
mentary operator Q = 1 — V . Then applying the iden- 
tities VCoe = £mP , VCqi = CoiV = we obtain the 
following set of differential equations for the functions 
u(t) = Vp(t) and w(t) = Qp(t) 

ii(t) = (£ 0E + V{Ci + C 2 )) u + V(d + £ 2 )w (D9) 

w(t) = Q(d + £ 2 )w + (Coe + Coi + Q(£i + £ 2 )) u, 

(D10) 



we obtain 



We solve (|D9[) - ([D10[) for u(t) using the Laplace transfor- 
mation. This leads us to 



u(t) = (C aE +V(Ci+C 2 ))u(t) 



r(d+C 2 ) dr 
Jo 



i-(C e+Coi + Q(Ci+£ 2 )) 



x Q(d + C 2 )u(t - t) 



(Dll) 

In the next step we apply the Markov approximation to 
the function under integral: u(t — r) sa u(t), and we 
neglect in (ID11I) all the terms of order higher than rj 2 . In 
this way we arrive at 



i(t) = (£ v+V(£i+£ 2 ))u(t) 



VCi dr 
Jo 



= t(£oe+£oi) 



QCiu(t) 



(D12) 



Finally we trace over the internal degrees of freedom ob- 
taining the ME for the reduced density operator p e = 
Tr/{p}. A calculation of the trace of the first term on 
r.h.s. of (ID12I) gives 



Tr,r{£oEw(i)} = CoEPe 
Ty^VCiu^)} = i^(a y ) s [zi,p e ], 
k 2 £l 

Tv^VC^t)} = [a x ) s [z 2 ,p e ] 

7,2, 



(D13) 
(D14) 



+ -k 2 (a + <r-) s (2zip e Zi - z 2 p e - p e z 2 ) , 
5 

(D15) 



where (.) s = Tij{.p s } denotes an average taken in the 
stationary state of the density matrix p s for the internal 
degrees of freedom. The calculation of the trace of the 
second term in Eq. (|D12[) is more involved, and as a result 



/>oo 

TniVd dre T ( £oE+£oi )Q£iu(t)} 
Jo 



k 2 n 2 r°°, 

dr 



{Zi{r)p e z l - ZiZi(T)p e ) 



x (<*y(T>„(0)) a -<<7„)2) +h.c. 



(D16) 



In the derivation of (ID16[) we have applied the following 
approximation e TCoE (zip e (t)) = (e TCoE Zi)(e TCoE p e (t)) ps 
(e TCoE Zi)p e (t), which is compatible with the Markov ap- 
proximation used in the above equations. The symbol 
Zj(t) denotes the operator evolved according to £qe' 



z 1 {t) 



e TCoE z- = e - ltH / h z . e itH / h 



while 



(a y (r)a y (0)) s = Tr z {a y (e TC °< a yPss )} 



(D17) 



(D18) 



is a two-time correlation function, discussed in more de- 
tails in the next section. 



Appendix E: Dynamics of the internal degrees of 
freedom 



In this section we analyze the dynamics of the internal 
degrees of freedom, which to lowest order in the Lamb- 
Dicke parameter is governed by 

p(t) = £oiP, (El) 

with £oi defined in (|D2[) . We start by writing Bloch 
equations for the averages of the ion internal operators: 



>xh ( a y)i an d ( a z), which follow from the ME (IE1I) 



dV 1 



A(cr T 



dt 



(E2) 
(E3) 
(E4) 



From Eqs. (|E2|) - (|E4|) one can easily calculate the average 
values of the internal operators in the stationary state 

(E5) 



S 2 ' 

7 2 + 4A 2 
4S 2 ' 



(E6) 
(E7) 



where 5 2 = ^ft 2 + \^ 2 + A 2 . Now we utilize the quan- 
tum regression theorem (see for instance Ref. |33|), 
and with the help of Eqs. (|E2|) - (|E4[) we write the equa- 
tions of motion for the two-time correlation functions 
s k (t) = (a v (0)cr fc (r )) , where k = x,y,z, 



s x (t) = As y - 



-fis 2 



As x 

- is z - 



7 



(E8) 

(E9) 
(E10) 



2 " 

Applying the Laplace transformation we convert the set 
of differential equations into a set of algebraic equations 
and solving for s y (s) = J °° dr e~ ST s y (r) we find 
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Sy(s) 



(s + 7 )(a + |K(0) - (s + j)As x (0) ~ 0(f + |) MO) - »7 2 /(2.5 2 S )) 
(s + 7 ) (A 2 + (s + ^) 2 ) +ft2 (s + 2) 



(Ell) 



The values of the correlation function Sk(r) at t = 
can be easily expressed in terms of the averages (|E5[) - 
& s a (0) '= -i<<T,)„ « y (0) = 1, «,(0) = iK) s . Now 
observing that 



dTe^ T (a y (T)a y (Q)) s = s y (iu>)*, (E12) 
we can calculate the function S(ui) defined in (|47|) : 



(E13) 

where 5(uj) is the delta function and Pf(l/a;) denotes the 
principal value. 



Here, u>i(t) — + 2Acos(u>ft)), 



w k 


= Im 


" poo 

/ dr fk(T)(<J y (T 
Jo 


^(0)} s 


v k 


= Rc 


" poo 

/ dr fk(r)((T y (T] 
Jo 


°v(0)>. 


u k 


-r 

Jo 


dr f k (r), 





(Fl) 
(F2) 
(F3) 



Appendix F: Calculations beyond RWA 



In this section we derive equations of motion for aver- 
ages quadratic in the operators Zi , pi , z m , p m . Assuming 
that r)Q <C uJi, we neglect the second and the third term 
in the commutator in describing the effects of the 

laser field on the ion's trapping potential. Introducing 
the dimensionless variables x\ — (z 2 )/l 2 , x 2 = { z m ) l^mi 

^3 = (P 2 )l 2 / h2 i X i = (P 2 m) l rn/ h2 i x 5 = (*iPi + P l Z l )/H, 
%6 = (z m Pm + PmZ m ) / 'K X 7 = (z m Pi}k / '{Kim) , ^8 = 
(zip m )l m /(Hli),X 9 = (ptPm) klm / K 2 , X 10 = (ziZ m ) / '(Mm) , 

we obtain the following set of linear differential equations 



x\ =ujiX 5 
x 2 =u) m x 6 , 

&3 



= - UJ l (t)x 5 + jJlf((T + (T-) ss - Auj c X 7 

- Wix b - 2W 2 x 3 - 2WiX 9 - 2W 3 x 7 - 



±4,= — LJ m X e - 4:UJ C X 8 

x 5 ~2bJiX 3 — 2uj i (t)x 
-W 2 x 5 



l 



Xq =2uj m X4 : — 2u} m X2 

x 7 =uj m x g - u}i(t)x w 



- 2W x Xi - 2W 4 x 8 - 2W 3 x w 

- 4w c xi , 

— 2u) c X2 

- W 2 x 7 - W3X2 - Vi 



-V2 + U2, 



±8 =^iXg - ii) m X\0 - 2uj c Xi, 
m 



-Ui, 



- x&) 



=UJiXg — UJ m X W - lUJ c Xi, 
X 9 = - UJ m X 7 - UJi(t)x S - oj c (x 5 

- \W 3 x & - W lXs - W 2 x 9 - W4X4 + V 3 - U 3 , 

X 1Q =L>iX 7 +UJ m X$. 



with k = 1,2,3,4. The functions fk{r) are determined 
by the evolution of the operator Zi due to £qe: Zi(r)/li = 



0TC0 



■'Ik 



fi(r)zi/k + f2(r)pih/h + f 3 (r)z m /l m + 



fi{r)p m l m /h. We note that Zi{t) = Zi(—t)u, where Zi(t)n 
denotes the operator evolved in the Heisenberg picture. 
The functions fk(r) can be easily calculated from the dif- 
ferential equations derived from £oe (|D3[) , with H given 
by (ET 



dpi , x 

— — =Uli(t)Zi + 2to c z m , 
dr 



dp r < 



dr 

dzj 

dr 

d7 

~dT 



= - UiPi, 



(F4) 
(F5) 
(F6) 
(F7) 
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